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Abstract. The recent anomalous segregation experiment ||l|| shows dramatic, rapid internal 
state segregation for two hyperfine levels of *^Rb. We simulate an effective one dimensional 
model of the system for experimental parameters and find reasonable agreement with the data. 
The Ramsey frequency is found to be insensitive to the decoherence of the supeiposition, and is 
only equivalent to the interaction energy shift for a pure superposition. A Quantum Boltzmann 
equation describing collisions is derived using Quantum Kinetic Theory, taking into account 
the different scattering lengths of the internal states. As spin-wave experiments are likely to 
be attempted at lower temperatures we examine the effect of degeneracy on decoherence by 
considering the recent experiment where degeneracy is around 10%. We also find that 
the segregation effect is only possible when transport terms are included in the equations of 
motion, and that the interactions only directly alter the momentum distributions of the states. 
The segregation or spin wave effect is thus entirely due to coherent atomic motion as foreseen 
in|§]. 



1. Introduction 

A topic of great interest in many particle quantum mechanics is the effect of coherence in 
mesoscopic systems. The recent experiment at JILA [|l]] shows transient highly non-classical 
behaviour, characteristic of such phenomena. An ultra-cold non-condensed gas of *^Rb is 
harmonically trapped and prepared in a superposition of two hyperfine states by a two photon 
7r/2 pulse. The superposition is allowed to evolve for a range of different times and the 
densities of the two levels measured across the trap; the data are then collected via absorption 
imaging. In the cold collision regime of the experiment the two states interact via slightly 
different S-wave scattering lengths; moreover, the states experience different Zeeman shifts 
in the presence of a magnetic bias field, and the combination of these effects produces an 
effective potential between the two states This differential potential was characterized 
using the Ramsey spectroscopy technique, and the segregation effect was explored for 
different potentials and atom densities. When the differential potential is constant across 
the cloud the motion is not observed. When the differential potential has a gradient across 
the cloud the atoms appear to redistribute in the trap, as long as the two states are in a partial 
superposition. The dissipation of the coherent relationship between the states is responsible 
for the transience of the spin waves. 

The original experiment stimulated much work on this system. Three theoretical 
treatments |[ ^] found good agreement between simulations of the segregation dynamics 
and the evolution of the experimentally measured distributions reported in [p. These works 
all use the spin operator notation and treat the system as essentially a spin wave problem. 
More recently, detailed imaging of the Bloch vector was carried out during spin wave motion 
in [M]. Simulations of linearized spin kinetic equations also showed good agreement for the 
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frequencies and damping damping times of linearized collective modes We do not use a 
spin operator formalism in this paper, although this has no effect on our results. 

The outline of this paper is as follows. The equations of motion are found from the 
Hamiltonian for the system in Section |^. The spectroscopy of the initial state is outlined 
in Section ||. The interpretation of the Ramsey experiment and the interaction energy are 
discussed in Section ^ In Section |] the issue of damping is addressed. A Quantum Boltzmann 
equation describing the effects of collisions on the distributions for states with different 
scattering lengths is derived, and we discuss the relaxation time approximations used in 
the simulations presented in Section ^. In Section ^ we consider the question posed in the 
experimental report regarding the physical state of motion of the atoms during the segregation. 



2. Theoretical Framework 

2.1. Second Quantized Hamiltonian 

The second quantized Hamiltonian for a coherently coupled two state gas is 
H= J cPx |7/4i7i(x)V'i + ^|i?i(x)V2 



+ 9E{^) (^Vl^2e*''*+V'|V'ie-*''*) 



where 



(1) 



Hj (x) = - ^ + V, (x) + hiu, (x) (2) 

for j = 1,2, and g is only nonzero during the initial two photon pulse, creating a coherent 
superposition of the two hyperfine states. Note that the Zeeman splitting of the transition 
frequency is position dependent but this is absorbed into the effective external potentials. The 
Heisenberg equations of motion for the field operators are 

ihi^i = ffi(x)V'i + uh^Z-ItAi^i + wi2V'2V'2V'i + .9£^(x)V'2e*"*, (3) 

ihlp2 = H2{-X.)lp2 + U22lpllJ-'2lp2 + Ul2'0l'0lV'2 + .9£'(x)'i/'l 6^*"* . (4) 



2.2. Experimental Details 

The experiment uses an ultra cold non-condensed *^Rb gas held in a harmonic trap at 
temperatures T~ 850 nK. The trap is a cigar shape with frequencies ujr = 27r x 7 Hz 
and LOz = In X 230 Hz. A 7r/2 two photon pulse is used to create a superposition of the 
|1) = \F — l,mp = —1) and |2) = \F = 2,niF = 1) hyperfine states from a thermal 
equilibrium |1) gas. The s-wave scattering lengths between the hyperfine states |i) and \j), 
aij, are an = 100. 9ao, ai2 = 98.2ao and 022 = 95.6ao, where gq is the Bohr radius [|l]]. In 
the usual cold collision pseudo-potential approximation this leads to the interaction strengths 

u,j = ATih^aij/mR. (5) 

where m is the *^Rb mass. Since the scattering lengths are so similar, and almost evenly 
spUt about we define the interaction splitting 

6u — uii — U22 = 5.3 , (6) 

m 
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and use the approximations 



"11 - "12 = 2.7 ~ du/2, (7) 

m 

ATTh?ar) , , 

"12 - "22 = 2.6 ~ 5u 2, (8) 

m 

in the remainder of this paper 

A magnetic bias field is used to change the effective differential force exerted on the two 
states, and can be chosen to either cancel or enhance the mean field force. Segregation of the 
two species is observed via a subsequent two photon pulse which causes a transition to an 
internal state appropriate for absorption imaging. The magnetic bias field is used to control 
the onset of segregation for a given density, or alternatively the density may be increased for 
a fixed bias strength to produce a similar effect. 



2.3. Equations of Motion 

In order to simulate the full behaviour of this system we wish to find suitable equations 
of motion to describe the system in the Hartree-Fock regime, where a local density 
approximation is valid. Defining the Wigner amplitudes 

n,(x,p)- yrf3y(V'](x-y/2)V,(x + y/2))e-'P•y/^ (9) 
/(x,p) = / o^3y(^'^(x-y/2)^'2(x + y/2))e-*P■''/^ (10) 



so that the densities are 

= /,/,t/' 



NA^) = J ".(x,p) ^ (V'](x)^,(x)), 



(11) 



-Vydiff(x) • Vpm(x,p) 



^(x) = J /(x, p) ^ (^I(x)V'2(x)), (12) 

and using Hartree-Fock factorisation and standard Wigner function methods the equations 
of motion are written in terms of rt(x, p) = ni(x, p) + 712 (x,p) and the segregation 
m(x, p) = ni(x, p) - (x, p) as 

n(x, p) = {-(p/m) • Vx + VK=ff (x) • Vp} n(x, p) 

_ 1, 
2 

+ "i2{VF(x) • Vp/(x,p)* + VF(x)* • Vp/(x, p)}, (13) 
m(x, p) = {-(p/m) • Vx + Vycff(x) • Vp} m(x, p) 

- ivydiff(x) • Vpn(x,p) 

- ^{/(x,p)^^(x)* - /(x,p)*F(x)}, (14) 

/(x, p) = {-(p/m) • Vx + Vyeff(x) • Vp} /(x, p) 

.Vdiff(x) 

- « — /(x,p) 

- (m(x, p)F(x) - /(x, p)Af(x)) 

+ ^Vf^(x). Vpn(x,p), (15) 



Theory of the Ramsey spectroscopy and anomalous segregation in ultra-cold rubidium 



4 



where the potentials are 

K=ff(x) = i(^V^i(x) + y2(x)+3ui2A^(x)+5uA/(x)^, (16) 

Vkiff (x) = ;i A (x) + HAlUc (x) . (17) 
and the frequencies are the coherent frequency shift 

/lAwe(x) = 2w22A^2(x) - 2wnA^i(x) + 2wi2(A^i(x) - iV2(x)) 

~ -faiV(x), (18) 

and the Zeeman shift 

;iAw,(x) = n(tj2(x) -cc;i(x)) = F2(x) - Vi(x). (19) 

For the low densities employed these equations of motion are nearly exact for the thermal gas, 
apart from the relaxation caused by collisions. 

3. Spectroscopy of the initial state 

The result of the laser excitation is to rotate the spin wavefunction into the |2) subspace, so 
the one-body wavefunction is transformed to 

|1) ^i|l)+r|2), (20) 

with |tp + |rp = 1. In field theoretic language (in the Heisenberg picture) we describe the 
initial state by the field operators ?Aji(x) and ipi2{x.), where the subscript i denotes refers to 
the initial state before the pulse. These transform to 

ipiiix.) -01 (x) = iV'ii(x) + rV'i2(x) (21) 

0i2(x) -02 (x) = -ripii{-K) + t'ljji2{x.), (22) 

where the condition |ip + jj'P = 1 ensures the conservation of total number and r and t 
are chosen real. We denote the population densities by Nj{x.) ~ (?/;j(x)-0j(x)), and the 

coherence amplitude by i^(x) — {ip\{x.)ip2{^))- In general the population densities transform 
to 

iV,i(x) ^ Ni{^) = |tpiV,i(x) + \r\^N,2{^)+t*rF,{x.)+r*tF*{^), (23) 

iVi2(x) -> 7V2(x) = |r|2iV,i(x) + \tfNa{^) - rri^*(x) - r*tF,(x), (24) 

For the ^^Rb experiment the initial occupation of the 2 state is zero, so that Nn (x) — N{x.), 
Ni2 (x) ~ Fi (x) — and a two photon pulse causes the transformation 

A^a(x) ^iVi(x)= |tpAr(x), (25) 

A^,2(x) ^iV2(x)= |r|2Ar(x), (26) 

F,;(x) ^ ^^(x) = -rrA^(x), (27) 

|F,(x)|^ ^ \F{^)\^ = 7Vi(x)7V2(x), (28) 

so that (Est) corresponds to a fully coherent superposition. 



3.1. Energy Density 

If we consider the total energy density we have in general 

t/(x) . U + V^I(X)) -0,) + (02 (-^ + y2(x) ) V.2 



(-y-'/'l^lV'lV'l + Ml2'0lV'20102 + ^V'2V'2V'2'02) • (29) 
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Changing the internal state of an atom does not change its velocity so the kinetic energy 
terms may be ignored in finding the change in energy of the gas during the transition. Using 
Hartree-Fock factorization and (|2^-p8|), the energy density is 

C/(x) = unNii^f + wi2(|F(x)|2 + iVi(x)iV2(x)) + U22N2{^)^ 

+ l/i(x)iVi(x) + y2(x)A^2(x). (30) 

Evaluating the change in energy density for an arbitrary rotation defined by (pH), (p2]) leads to 



A;7(x) = ^unilt - 1) + 2^12^1 VI' + U22\r\^jNi^f 

+ iV2-Vi)\rfNi^). (31) 

3.1.1. Infinitesimal rotation For a very small rotation |rp ^ 1 and 1 — ~ 2|rp. Putting 
A^2(x) = SN{x.) < A^(x), |r|4 ~ and iVi(x) ~ A^(x), we find 



A(7(x) = SNi^) (^2(ui2 - Un)iV(x) + V^2(x) - V^i(x) ) , (32) 
For ^^Rb this becomes (using (^ ^) 

A(7(x) = SN{^) ( - Su 7V(x) + y2(x) - Fi(x) ) = (5Ar(x)Kiiff (x), (33) 



so that the change in energy density varies linearly with the number of atoms transferred to 
the higher energy state, and is proportional to the differential potential. 

3.1.2. tt/2 pulse For a 7r/2 pulse the energy density change is 

A[/(x) = (^3un + 2^.12 + ^^22) + (V^2(x) - V^i(x))^. (34) 

For *^Rb this becomes 

A;7(x) = (sim + Su) + (y2(x) - Fi(x))^. (35) 



Thus the effect of an intense pulse is to produce a nonlinear change in the energy density via 
the mean field interactions. 

3.2. Coherence energy 

The analysis is simplified by noting that because the only net force in the experiment is a weak 
differential force between the two internal states, the total atomic density is approximately 
conserved. The experimental data suggests that this is a reasonable approximation [|l]], and 
our detailed simulations confirm this expectation. Treating A^(x) as constant allows insight 
into the energy conserving processes in the dynamics. In terms of the normalized moments 
(p^(x)) = / dPp p''nj{x,p)/Nj{x){27rh)'^ we have the relations 

7Vi(x)(pi(x)) +7V2(x)(p2(x)) = 7V(x)(p(x)) = (36) 

(which is essentially momentum conservation in the absence of any external perturbation of 
the stationary density profile), and for the total kinetic energy density 

iVi(x)(pf(x)) ^ jV2(x)(pi(x)) _ jV(x)(p^(x)) ^ ^^^^ 
2m 2m 2m 
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Here only N{x.) is constant since the other moments and densities will change during the 
relative motion of the two species. We may now write in the form 

C/(x) . '-^Ni.f + I^^M±^A.(x) (38) 

+ ^M(x)2 + ^7V(x)A/(x) - i?iAc^,(x)A/(x) (39) 

+ wi2|F(x)|2 (40) 

+ iV(x)iP!M, (41) 
2m 

where the spatial segregation is = J (Pp m(x, p)/(27r?i)'^. The first line is constant 

in time, and the second line only varies with Af(x). The third line is the coherence energy 
density which may change during the motion, while the last line is a function of the local 
densities and temperatures through (|37|). It is clear that the coherence energy moves between 
F(x) and A/(x) via relative motion which changes the local temperature and momentum of 
each distribution. 



4. The Ramsey frequency and the interaction energy shift 

The Ramsey technique has recently been shown to be a particularly useful tool for exploring 
the coherence properties of ultracold dilute gases [|H pTj ^J. In applying the technique 
to the ultra-cold non-condensed gases of these experiments there remains an issue in the 
interpretation of the technique which we wish to resolve. The problem is whether or not 
the Ramsey frequency is sensitive to the decoherence of the superposition of internal states 
required to resolve the Ramsey fringes. The Ramsey technique has been developed and used 
extensively for probing the relative energies of internal states in non-interacting atomic beams. 
When interactions are negligible the measured quantity is the energy difference between the 
two internal states of the superposition. We will see that for the interacting case, this is only 
partially true because interactions produce a transient mean field energy which decays with 
the coherence of the superposition. 



4.1. Experiment 

The experiment measures the frequency of oscillation for the relative occupation of the two 
internal states, after a second tt/2 pulse is applied to the system. The final two photon pulse 
transforms the fields to 

ipi -> — (42) 

; i^^-i^i ,,,, 

-> (43) 

In terms of the densities Ni{x.) = and the coherence ampUtude F(x) = ('(/'i'02), the 

ratio of the densities may be written 

7V2(x) ^ jV(x) - 2Rc f (x) 
iVi(x) iV(x) + 2Re i^(x) ' 

where A^(x) = ^i(x) + A^2(x) may be assumed invariant. The ratio will then only depend 
on -F'(x); in particular, it will oscillate at the same frequency. Applying the second 7r/2 pulse 
at different times and measuring A^i(x)/iV2(x) thus determines the frequency of -F(x), and 
this is the frequency accessible via the Ramsey technique of | 
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4.2. Coherence equation 

In a typical experiment ^ ||], a tt/2 pulse generates an equal superposition of the two 
hyperfine states which is then allowed to evolve without coupling until a second tt/2 pulse 
is applied. We are primarily interested in the influence of the interactions on the Ramsey 
frequency, and after the first pulse the Heisenberg equations of motion generated by the 
interaction terms in the Hamiltonian ([l|) are 

ihlpi = Uiil/'iV'l'/'l + Ul2lp2'^2'4'l, (45) 
ihl]j2 = U22i'2i'2'^2 + Ml2'0l'0lV'2- (46) 

The operators t/'Jt/'i and ili\i)2 commute with the Hamiltonian when i?(x) = 0, so the 
occupations are preserved. The equation of motion for the coherence is 

'^^^f^ = U22{4'14^2^2) + Ui2{4'14M2) 

- Wll(V'IV'lV'1^2) - Ui2{lpltpltp2i^2)- (47) 

4.3. Four point averages 

There are two equivalent ways of treating the four point averages when the system is 
noncondensed and thermal, both of which lead to results of the form 

=2iVi(x)F(x). (48) 

i) When the initial quantum state is nondegenerate, each field may be written in terms of a 
set of orthonormal single particle wavefunctions 0^ (x) 

ipi =^(t)r{:>i-)air, (49) 

r 

where [a^r, cl^ji^] = SijSrk- The densities and coherence then read 

{4^^,) = E M^rM^){aU,s). (50) 

r.s 

The eigenvalues of the operators Nir = al^air which contribute significantly are either 
or 1 . We can describe the situation by the averages 

{ala,r) =N„, (51) 

{a\^a2r) = Mr. (52) 

For all averages {a\^a.js) = when r ^ s, corresponding to independently occupied 
modes. We then have 

r 

and the four point averages become, for example 

, t t X I SrjSskNirMs + 5sj5rkNisMr when r 7^ s and j 7^ fc, 

{a[ra[saija2k) = \ ^ ■ (54) 

I otherwise. 

The four point average becomes 

{44^2)^2 \M^)\^\M^)\^NirMs, 

r,s. r^s 

~27Vi(x)F(x), (55) 
where the approximation is valid if the occupation is spread over very many modes, as 
for a thermal state, so that the term with r = s is negligible. 
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ii) Alternatively, the same result may be found using Hartree-Fock factorization for the 
averages over the field operators. This method uses the Gaussian statistics of the thermal 
gas, which may hold under more general circumstances than the arguments of i), but is 
equivalent for the system under consideration. The averages are 

{^I4t^3^4) = {^%){4i^4) + {^\^4){4^3) + {^14) (^3^4)- (56) 

For a thermal gas there is no anomalous average, so we recover the same result as ( ^5| ) 

=2iVi(x)F(x). (57) 

4.4. Ramsey frequency 

Using either approach, the equation of motion for the coherence becomes 

ifi^F(x) = [2u22A^2(x) - 2uiiiVi(x) + 2ui2(iVi(x) - A^2(x))] ^^(x) 

= niUR{^)F{^). (58) 

Thus if the gas is thermal, the frequency of oscillation lor{x.) is independent of the amplitude 
|F(x)|. In particular if some kind of damping reduces |i^(x)| with time, this does nof change 
the Ramsey frequency. 

4.4.1. Transport and trap effects The full equation of motion (|5]) may be integrated to find 

+ V • (F(x)vj^(x)) = - I (Ac^e(x) + Aa;,(x)) F(x) 

at 

= - ^u;i^(x)F(x) (59) 

where the velocity is 

1 r d^p 

toF(x) J {27rhf 

When the velocity vanishes we recover ( p^ with an additional shift caused by the trap energy 
difference for the two states. The coherence current described hy vp may alter the measured 
Ramsey fringes, and indeed full simulations may be required for comparison with experiment 
when significant motion occurs. 

4.5. Interaction energy shift 

The quantity which is sensitive to the loss of coherence is the interaction energy, which in 
general does not correspond to the Ramsey frequency Wfl (x). 

The interaction energy change caused by changing the internal state from |1) to |2) is 
found from the chemical potentials associated with this change 

Neglecting the transport and trap terms for simplicity, we use the mean energy density for the 
interactions 

C/(X) = +u,2{44M2) + ^(44^24^2), (62) 



v^^W = I 7:^P/(^'P)- (60) 



and factorize averages as in (55 1, to find 



g\p\2 

Hi = 2unNi + U12N2 + (63) 

oNi 

Q\p\2 

(12 = 2U22N2 + U12N1 + -j-^. (64) 

0N2 
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The change in energy caused by this transition is A/x — — l^i = ^Aw^, so that using (|63[), 
(g) we find 

ftAcj^ = 2U22N2 - 2uiiNi + ui2{Ni - N2) 

The derivatives of the coherence energy density are not evaluated because it is not usually 
possible to do so in any direct manner since |F(x)p cannot generally be specified by 
knowledge of the -/Vi(x). The Cauchy-Schwartz inequality leads to 

|F(x)|' < 7Vi(x)Af2(x). (66) 

When |i^(x)p ~ Ni{x)N2{yi) the superposition is purely coherent, leading to a factor of 2 in 
the cross interaction part of the frequency shift; whereas for \F\'^ = the extra factor due to 
the coherence is absent and we recover the thermal result. A simple model which we will use 
is found by taking |Fp ~ a{t)N\N2, with a(0) = 1 reconstructing the initial condition (^8|), 
and a{t ^ 00) = medelling the damping to thermal equilibrium. The mean field energy 
shift becomes 

ftAcj^ 2U22N2 - 2MiiiVi + 2a{t)ui2{Ni - N2) (67) 

Immediately after the pulse 

Atj^(x) = WflXx), (68) 

so that for short times the Ramsey frequency coincides with the interaction energy shift. 
Inclusion of the effect of collisions will cause |F(x)| to decay with time. The simplest way 
of doing this is to consider the case where there is no differential potential gradient so that 
segregation cannot occur The collisions will still have a relaxation effect on the coherence 
and as a first approximation we may use the equation of motion 

F = -iojR{x)F - —F, (69) 

Tr 



where Tc is the timescale of relaxation. Solving this and using (65 ) gives 



hAuj^ ^ 2U22N2 - 2unNi + ui2{Ni - N2){1 + c-^*^^^) 



(70) 



Clearly the mean field energy shift is sensitive to the decay of coherence, and decays at twice 
the rate for the amplitude. 



5. Kinetic Theory of Coherence Damping 

We now consider the role of collisions in dissipating coherence. Our starting point is the 
quantum kinetic theory of Gardiner and ZoUer [^], which may be used to find an expression 
for the damping rate of /(x, p). While the expressions found are quite explicit, they are not 
easy to simulate, so that for the actual simulations we will use simplified forms based on a 
relaxation time approximation. 



5.7. Quantum Kinetic Theory with internal degrees of freedom 

In order to derive the damping for /(x, p) we first consider the damping term in the equation 
of motion for the reduced density matrix for a gas with a single internal state. For brevity 
of notation we suppress x arguments where possible since the damping is local. For ease 
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of comparison with [^] we carry out the following in K variables. We seek the equation of 
motion for the terms 

/,,„(x, k) = / dV(^'J(x - y/2)^„,(x + y/2))c-^''-y (71) 



where l,m = 1,2. The derivation of the Quantum Boltzmann equation may then be carried 
out along similar lines to that in [^. The principal modification of the derivation is due to 
the different scattering allowed between distinguishable internal states of the particles, and 
the preservation of possible coherence between internal states. The resolution of the field 
operator for internal state i into its momentum components V^jk (x) is defined similarly to 
(31) of 01 as 

V'^Cx) =^e-*K.x^,K(x). (72) 

K 

The projector defined in (48) of for the derivation of the master equation simplifies slightly 
because no Bose-Einstein condensate is present in this case, and becomes 

I otherwise. 

This projector identifies all configurations with the same distribution in K, but leaves the 
position and internal state distributions undisturbed. To evaluate the effect of collisions on the 
coherence we require the equation of motion for the reduced density matrix 

"Pn/O = PnPPn = wn, (74) 
induced by the interaction Hamiltonian 

^/ = ^E^'j(l'2.3'4), (75) 

ij, e 

in which we have used the notation for the sum over all momenta defined in [0] 

e Ki,K2,K3,K4 

and the interaction operator [/^ (1234) = C/y (e): 

[/y (1234)= j d^x j d^x' c(»Ki x+,:K2 x'-»K3 X-iK4 x') 

X V'Iki WV'lKj^O'^ylx - x')^»K3(x')V'jK4(x). (77) 

In this form Hj is just the K space representation of the second line of (|I]). In what follows 
we will make the usual psuedopotential approximation for the interaction parameters and set 

(x - x') = '^(x - x') = Uy (5(x - x'). (78) 

In order to describe the exchange collision 1, 2 ^ 4, 3 we also define the interaction operator 
C/ij(1243) = ?7y (e) as 

U,^(e) = U., j d-^X c'(Ki+K.-K3-K4).x 

X 4k, (x)^']k. (x)V',K3 (x)V',;k4 (x) (79) 
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which reverses the momenta of the last two operators. Carrying out the procedures of 
leads to the density matrix collision term similar to (75e) of 



coll 



^ E ^(A-(e)) 



ijlJ. e 



+ 



U,j{e)vN-oit)Ujj{e) + Uij{e)vN-eit)uUe) 



(80) 



where Au!{e) = uj^ + uj-^ — u!2 — The first line is the contribution from the collisions 
1,2^3,4 and 2, 1 ^ 4, 3 and is a sum of terms like (75e) for a single internal state , while 
the remaining terms describe the collisions 1, 2 ^ 4, 3 and 2,1-^3, 4. Of course if only one 
internal state occurs, the substitution i — j = I = J restores the exchange symmetry of the 
operator Uii{e) = Uuie) = C//j(— e) and we recover the ordinary collision term. 

5.2. Modified Quantum Boltzmann equation 

We now seek the Quantum Boltzmann equation for this system, which is the equation of 
motion for (71 ) in the continuum limit. Tracing over the reduced density matrix allows us to 
write the equation of motion for 



(V'/k'^'nk) = y^Tr{t!NV'/kV'mk} 



(81) 



N 



in the form 



ijlj. e 



X I ((7|,,(e)[7^/k^,„k, U,, (e)]) + ([C/y (e), ^./kVAnk]C/|,,(e)) 



(82) 



+ (;7|,,(e)[<V^™k,C/.,(e)]) + (e),<V™k]C/I,(e)) 



(83) 



The terms in lines ( |82| ) and ( |83[ ) give the direct and exchange scattering contributions 
respectively. The details of passing to the continuum limit are discussed in A^ppendix A. 
For the purposes of ease of comparison with other calculations of this kind and coherence 
with the rest of this paper, we will also write the final result in terms of momenta P = KK 
and energies e = huj instead of frequencies and wave vectors. We eventually find 



dt 



coll h J {2TThf J {2TThf 

X 5{P + P2 - P3 - Pi)5{€ + £2 - £3 - £4) 



ilj 



UIJ 



ifuiP) + Slj){fu{P2) + 5u)fUP3)f.HiP4) 

- /ij(P)/./(P2)(/7m(P3) + (5/™)(/j»(P4) + Sj,) 
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+ (//j(P) + 5u){fu{-P2) + (5,7)/7,(P3)/j™(P4) 

- /zj(P)/./(P2)(/,7™(P3) + <5j™)(//,;(P4) + 5n)^u, 
(^(/jm(P) + <5j™)(//.(P2) + <5/0/;/(P3)/a(P4) 

- /jrn(P)/7.(P2)(/i/(P3) + <5,7 ) (/,; J (P4) + <5a) 

+ (/j™(P) + <5j™)(//,(P2) + <57.)/«/(P3)/ij(P4) 

- /j™(P)/7^(P2)(/l/(P3) + 4/)(/./(P4) + -5,/) 



(84) 



This is a sum of direct and exchange Quantum Boltzmann collision terms, summed over 
the possible interactions for each case. This is now in a form where simplifications may be 
introduced, and there are a number of useful results that can be obtained from this equation 
for the particular case of ®^Rb. 

5.3. The case of equal scattering lengths and low density 

If we drop the terms proportional to 5ij, and set the scattering lengths equal, so that Uij = ui2, 
we straightforwardly recover the Boltzmann limit for / 7^ m, in agreement with (6) of [^, 



dt 



h J {2TThf J (27rft)3 
X 6{P + P2 - P3 - P4)'5(e + £2 - £3 - £4) 

X {3n(P3)//™(P4) + 7l(P4)//™(P3) - ?l(P)/i™(P2) - 3n(P2)/,™(P)} , 

(85) 

where we have set /n + /22 = n. These terms involve products of two distributions, and 
we will use a notation in which terms involving products of j distributions will be denoted by 

flm \j ■ 

It is important to note that the presence of coherences in the collision integrals means 
that the usual arguments for the relaxation time approximation are not strictly applicable. For 
the 1 1) state for example, with the same approximations as above, the collision term becomes 



df 



11 



dt 



d^Po f d^P 



d^p. 



h J {2TThf J {2TThf 
X <5(P + P2 - P3 - Pi)5{e + £2 - £3 - £4) 
X {2n(P3)/ii(P4)-2n(P2)/ii(P) 
+ 2/n(P3)/ii(P4) - 2/n(P2)/ii(P) 

+ 2/l2(P3)/2l(P4) - /l2(P)/2l(P2) - /l2 (P2)/21 (P) } , (86) 

SO that we expect the coherence terms, like those in the last line, to provide extra damping for 
the distributions. 



Theory of the Ramsey spectroscopy and anomalous segregation in ultra-cold rubidium 13 



5.4. Equilibrium coherence damping 

The Quantum Boltzmann collision integrals for the distributions n; (x, p) vanish if the 
distributions are in local equilibrium. When there are two interacting internal states with 
different S-wave scattering lengths, there is an additional damping effect caused by the cross 
interactions. In the experiment the degeneracy is typically about 10% so we will find the 
effect of keeping the differences in scattering lengths i n (p4| ) using the Boltzmann equilibrium 
momentum distribution. Since the quartic terms in (^4[) cancel, there are two processes to 
consider. The quadratic terms provide the most important contribution for low phase space 
density. The cubic rate has the additional interesting feature of depending on the phase space 
density A^(x)Aj[j, where the thermal deBroglie wavelength is Ath = {27rh? /mkBTy/^ . 

5.4.1. Quadratic terms Retaining the differences in scattering lengths we find the term 
corresponding to the Boltzmann approximation (p5|) is 



Irn 



dt 



d^P4 



2 h J {inhf J (27rft)3 

X 5(P + P2 - P3 - P4)<5(e + 62 - £3 - £4) 

(/„n(P3)/i,(P4) + /«(P3)//m(P4)) 

- 4(/™.(P)/;,(P2) + flmiP)UP2)) 

- «L(/H(P)/^m(P2) + flm{P)UP2)) ^ (87) 

To find the local damping rate of the coherence amplitude A^/,„(x) = J ■j^^^^jafimi'^, p) we 
approximate the equilibrium distribution by the local Boltzmann equilibrium form 

(Q/ \ 3/2 
-J exp(-ap2), (88) 

with a = l/2mkBT. Integrating leads to 

dNi, 



dt 



hi^ II i'^a - u„nf (Ar,;,„(x)iV,,(x) + A^„;(x)Ar;™(x)), (89) 



where 134 is calculated in |Appendix~A Clearly the collisions between the same hyperfine 



state do not cause damping of coherence, since this expression vanishes for I = m. When 
/ ^ TO, the interaction terms are {uu — Ui„i)^ = {5u/2)'^ for all i. In terms of the thermal 
relative velocity, the final result may be written as 

9F(x 
dt 

where 5a = an — ai2- The rate coefficient takes the form avp{x.), which occurs in S- 
wave scattering of distinguishable particles, but with an effective scattering length Sa. For the 
experiments of [|l]] this produces a 1/e decay time for the coherence of ~ 13 s at the center of 
the trap, as noted in . 



= -ATT(SayVrN(x)F{:K.), (90) 

2 



5.4.2. Cubic damping For single species collisions, the cubic terms in the Quantum 
Boltzmann equation are responsible for the Bose enhancement of scattering which is crucial 
in describing BEC growth, four wave mixing and other statistical collision phenomena that 
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can occur in nonlinear atom optics at high phase space density. For two internal states, the 
cubic terms of (p4|) are 



dflr. 



dt 



3 h J {2T:hf J {2Trhf 

X S{P + P2-P3- P4)S{e + 62 - 63 - 64) ^ 

i J 

/ij(P)/^m(P3)/,/*(P4)«./ + h.j{P2)fjm{P3)lH{Pi)uu 

- f.Jm{P?^)fH{P)hj{P2W, - f.H{Pi)fu{P)lnn{P2)u^.J 

+ /^J(P)/^^(P3)/Jm(P4)^^J + kj {P 2) f .h{P 'i) fhn{P i)uij 

- fjm{P3)fu{'P)k(P2W, - f.H{Pi)fhn{P)kj{P2)u^.j\u, 



fjm{P)fH{P3)lu{Pi)uu + fM{P2)fu{P3)fun{Pi)u^U 

- flJ{P3)f^miP)fJ^iP2)uu - /. J (P4)/,/rn (P)/h (P2 )«i J 

+ fjrn{P)MP3)flj{P4)uu + /j. (P2)/,: J (Ps)/;™ (P4)^^™ J 

- fljiP3)fjmiP)MP2)uu - f.uiPi).fl,niP)fMP2)uij]uu 



(91) 



Using a thermal Boltzmann distribution and the same procedures as in the previous section, 
this may be reduced to 



ON,, 



dt 



J/234 



{Uiiuji + UimUjrn " 2Uim.U.Jl)iNiiNij Njm. + NijNjiNhn) 



(92) 



As for the quadratic rate, when I = m the cubic contribution vanishes. When I ^ mwe may 
use I23A from Appendix A to find 

9F(x) 



dt 



47r(H'^^.A?h (^i(x)2 + 7V2(x)2 + 2|F(x)|2) ^^(x), 



(93) 



where again we have separated the effective cross-section. For the rubidium experiment where 
degeneracy is about 10% we estimate the peak damping rate by simply using Nij ~ N/2 
which holds just after the first Tr/2 pulse. Equation (P3|) now becomes 



OF(x) 



dt 



' -47r{SafvrXl,N{x.)N{x.)F{:x.). 



(94) 



This rate is proportional to (Twp(x)A^[jp(x), which is distinguished from ( pO[ ) by an extra 
factor of the phase space density. For the JILA experiment [|l|] this process leads to a 1/e 
decay time at the trap center of ~ 240 s. The cubic rate is clearly unimportant in this regime, 
but if the phase space density is increased significantly, the cubic effects can approach or even 
exceed the quadratic damping. 
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It is also interesting to note the presence of a self interaction term in the cubic damping 
rate, proportional to |i^(x)p. It was found in [ pl| that the coherent velocity changing 
collisions responsible for spin waves may suppress the decoherence caused by atomic motion 
in a non-uniform system. For high phase space density there is also the possibility of the 
opposite effect becoming important, so that even in equilibrium the coherent interactions may 
cause significant additional damping. 

5.5. Damping to local equilibrium 

The equations of motion must be modified to include the effects of collisions on the 
distributions. This causes a relaxation of the momentum distributions towards equilibrium 
which will be dealt with using a relaxation time approximation. The rates derived from simple 
collision time considerations require modification to account for the effect of coherence on 
the relaxation. The coherence distribution between the two fields has a significant effect on 
the relaxation process because it plays a similar role to the single species distributions. 

5.5.7. Relaxation time approximation The collision damping term in the equations of motion 
for a single internal state is written in the relaxation time approximation as 

n(x,p)|,on = - r-i(x)(n(x,p) - n(^«'(x,p)) (95) 

where n'-^^^ (x, p) is the local equilibrium, and is damping time which may depend on 
position and the relative velocities of the two internal states. An estimate for the collision rate 
may be found from elementary kinetic theory if the correct quantum mechanical scattering 
cross sections are used. The hard sphere model of interactions leads to the collision frequency 
for a given particle of internal state i with particles of internal state j and density Nj (x) given 

by [| 

v,j = Nj (x)7r4 {vr{^)),, (96) 

where dij ~ {di + dj)/2 is the average hard sphere particle diameter and (iv(x))y is the 
mean relative velocity of particles i and j. If the masses are equal, that is mi = rrij = m, 
then in the absence of any net relative velocity for each state, the velocity assumes the thermal 
value 

Vr = ^JlQkBTlTTm. (97) 

The quantum formulation of the cross section is then given by the substitution ttc?^^ ~ dij, 
where try = 47ra^j (1 + Sij). It is convenient to define the local equilibrium distributions as 

n|7'(x,p) = 7V,(x) cxp (-(p - ((p,) + (p,))/2)VmfcB(T, (98) 

where Ti = Ti{x) and (p;) = (pi(x)) are average local temperatures and momenta. 
Summing over the contributions from both internal states gives 

n,:(x,p)|coH = - («,.(x)),;,a,;,(iV,(x)/2)(n,(x,p) - n|f '(x,p)) (99) 

- (w,(x)),,a„-A^j(x)K(x,p) - ^(x,p)) (100) 

where the factor of 1/2 is necessary to prevent over counting when determining the total 
collision rate between identical particles. This form of damping conserves the collision 
invariants and also allows a further simplification. From the definition of the moments we 
have the identities 

7Vi(x)(pi) + 7V2(x)(p2) = A^(x)(p) ~ 0, (101) 
iVi(x)(p?) + iV2(x)(p2) = 7V(x)(p2). (102) 
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so when the densities remain approximately similar it is then reasonable to set (pi) + (P2) = 
0, and make the replacement Ti + T2 = 2T in (|9^. If we assume that each momentum 
distribution is never far from its own local equilibrium we may also neglect (p^), so that the 
damping only arises from colUsions between distinct states. This leads to 

"T-(x, p)|co;/ = 

- (t;.(x)) 12^12 (A^2(ni(x,p) - n[l^\^,p)) - 7Vi(n2(x,p) - ni\'\^,p))^ (103) 

which may be written as 

M / ,u / , , ^ cxp(-pV2mfcBT) 
m(x,p)|coH = - (Wr(x))i2Cri2— — m(x, p) - Af (x) \ 

when the term proportional to n(x, p) — ri^j"*' (x, p) — ^2i^''(x, p), which is very small, is 
neglected. This means that the segregation momentum distribution relaxes at approximately 
the same rate at which each distribution would relax under the influence of collisions with 
the other internal state. Reducing this to an effective one dimensional rate requires that the 
distributions be radially averaged. The result is simply the above form divided by the circular 
transverse cross section T:{2^Jk^T/rnLo^Y . Using the thermal relative velocity, the colHsion 
term is now simply 

9 9 I m , I , ^ , X exp (— p^/2mfcBT) \ 
fn{x,p)Ui = -2a\2Ujl^-^N{x) \ m{x,p) - M{x) ^ ^^^J^^J - 1 • (105) 

The prefactor determines an effective segregation relaxation time at the center of the trap of 
^ 30ms for the experimental parameters. 

If there are coherences present between the two components the above arguments are 
not complete because the coherence amplitude provides an additional damping effect. It is 
not a priori clear what sort of relaxation approximation should be applied to the coherence, 
although the Cauchy-Schwartz inequality for the system leads to 

|F(x)|2 <7Vi(x)Ar2(x), (106) 

so that it appears reasonable to damp the coherence with the same relaxation time 
approximation used for the ordinary distributions. In our simulations we simply damp 
all /i7Ti(x, p) at the same rate, but we modify the rate to find the best agreement with 
the experiment, since at this level of approximation the damping rate is effectively a free 
parameter of the model. We find increasing the damping rate by a factor of 2 gives 
reasonable agreement with the experiment. Better agreement could be found by optimizing 
this correction for each experimental run, but this not been pursued here. 




5.5.2. Effect of relative velocity It is apparent from the experiment that in some regions of 
the gas the relative velocity of the two states becomes comparable to the thermal velocity 
during segregation. We treat the effect of relative velocity between the two distributions 
by including this into the calculation of the collision time t^. The average relative velocity 
between particles in state i with those in j at x is 

/ f W f ^3 f ,3 , »»(x, p)r^3^(x, pOIp-p'I 

(-.(x)>,- -J dp J dp ^^^(^)^^.(^) (107) 

Using Q = p + p', q = p — p'; assuming the momenta are distributed according to 

n,(x, p) ^ iV,(x) cxp {-a{p - p,(x))2) (a/nf^^ (108) 
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Figure 1. g(x) and its asymptote. 



where a = l/2mkBT, and choosing the mean local momenta to be aligned with the x axis 
Pi(x) = pi{x)'k, the average is 



8m 



■ exp {-aip,^ +P3^)) J d^Qexp {-aQ^/2 + aQ ■ {p, + p^)) 
J d^q |q| exp + aq • (p^ - pj)) , 



(109) 



where the temperatures of the two components are treated as equal. Since the integrals are 
spherically symmetric 

(vrix)).,, = ^^^e-"(P'+Pl)/Q(a/2,a|p, + p,|)/Q(a/2,a|p,-p,|), (110) 



where Iqia, h) and Iq(a, b) are discussed in Appendix B. In terms of the average thermal 
velocity (67]) the final result may then be written as 



where 



{Vr{,x))i 



l{x) 



2x 



x+(l + 2a;2)- 



-eri(x) 



(111) 



(112) 



Since g(0) = 1 this reproduces the standard result for a single component gas, and for the 
case where there is no mean relative motion of the two internal states. In the high relative 
velocity limit obtained by setting Vj(x) = — Vj(x) = xu(.x) with Vr ^ v{x), the large x 
limit g(a:;) y/ir/Ax yields 

{vr{x)),j^2v{x). (113) 

Figure |] shows the simple interpolation given by g . The effect of relative veloci ty on the 
relaxation may now be included by substituting (111) for the thermal result used in ( 104 ). 
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5.5.3. Variation with temperature Since g(2) ~ 2, when the axial relative velocity is 
\/TTVr ~ 1.1 Vr the damping rate will be doubled. For the JILA experiment the modification 
is somewhat less because the thermal velocity is of the order Vr ~ 20 mms~^, whereas the 
relative segregation velocity may be estimated to be at most ~ 15 mms~^, corresponding to a 
50% increase in damping rate. This is expected to be important for the spin wave experiment 
since coherence damping has significant effects on the motion, and the coherence damping is 
also modelled using a simple relaxation time approximation. Moreover, if the experiments are 
attempted at higher temperatures this effect becomes more significant because the timescale 
for initiation of segregation scales like ~ T^i/^ and the cloud length scales like ^ T^l"^ so 
that the segregation velocity must behave like ~ T^^^, whereas the thermal velocity scales like 
~ T^/^, so the ratio of the segregation velocity to the thermal velocity will be proportional to 
~ Clearly this correction will become relatively less important for lower temperatures. 

To find the effect of reducing the temperatur e on the coherence damping, we may 
comp are the temperature dependence of the rates (100) and ( po| ) with (^4|). In particular 



lOOj ) and ( |90| ) are proportional to ^kBTN{x.). The ratio of cubic to quadratic relaxation 
rates will vary as NT-^/"^, so that this process is expected to become more significant at 
lower temperatures. Since the peak density also varies as T^^/^ the ratio will vary as T^^. In 
particular, if the temperature is reduced by a factor of 3 for the experimental parameters, the 
quadratic and cubic rates ( pO[ ) and (|4|) will become similar The cubic collision terms (^Tj) will 
be important in this regime, and in the simplest approach, the relaxation time approximation 
(IOC ) would need to be modified to account for the relaxation caused by (^l|), with Uij ~ ui2. 



One could also carry out moment equation calculations similar to those of Nikuni [g] for the 
cubic contribution to the damping in linearized spin equations, but this not pursued here. 

6. Simulations 

To simulate the experiment the equations ( p3[p^ , p3| ) are reduced to an effective one 
dimensional description, and damping is included by adding the simple relaxation 



approximation term (105) to the resulting equations of motion. 



6.1. Simplifications 

The high anisotropy of the trap and low densities used allow the a reduction to a quasi- 
one dimensional model starting from a noninteracting initial condition. We may also take 
7V(x) = as seen in the experiment, and verified by our simulations. 

6.7.7. Initial motion The state of the gas before the 7r/2 pulse, given by the Boltzmann 
distribution 

n(x, p)|,.o = -^exp - ^ — j (1 14) 

for the single internal state 1 1), is the stationary solution of the equation of motion 

n(x, p) = {-{p/m) ■ Vx + V(Fi(x) + 2u^^N{^)) ■ Vp} n(x, p). (115) 

We may therefore neglect the density dependent effective potential 2itii A^(x) if 



2uiiViV(> 



Vli(x) 



2uiiiV(x) 2uii7V(x) 

^ ^ < , ^ ' < 1. (116) 



/ssT + 2uiiiV(x) - RbT 
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For the densities used in our simulations 2uiiN{0)/kT ~ 0.02 allowing the use of the 
noninteracting initial condition. After the pulse the system consists of an equal superposition 
of two internal states in identical external configurations, so we simply use 

,.(x,p)|«.Arcxp(-^-^), ,117) 

m(x,p)|t=o = 0, (118) 
/(x,p)|t=o = in(x,p)|t^o- (119) 

6.7.2. One dimensional equation The high aspect ratio of the trap (33 : 1) means that there 
is a separation of timescales, leading to a physical transverse averaging of the atomic motion. 
This means that we may average the distributions over radial coordinates and momenta 
and find a one dimensional equation of motion. The only change is that the interaction 
strengths are divided by the transverse cross section of the sample given by -kE?^ where 
the transverse radius R±^ is two standard deviations of the equilibrium density. A typical 
interaction dependent term has the form 

^^^|^ = ^.G(x)/(x,p). (120) 

All of the terms can be seen to be either of this type, or to have no interaction parameter and 
therefore to be unaltered by the radial integration. Putting 

/(x,p) = (^^) cxp {-muliy^ + z')/2kBT) 

f 1 



c^v{-{pl+pl)/2mkBT), 



\2TrmkBT ^ 

G(x) - Gix) (^J^ j cxp {-mcvUx'- + y^)/2mkBT), (121) 
and integrating over radial momenta and coordinates leads to 

Ot t:{2^ TTkBT/mulY 

where the denominator is clearly the transverse cross-section of the sample, at a radius of two 
standard deviations. The effective scattering parameter has the dimensions of energy x length 
required for a one dimensional interaction parameter. 

6.2. The differential potential 

The differential potential described by the parameter z/diff is defined by the experimenters [[l|] 
as 



^diff T-X — {—r^ ), (123) 



27r y m \ dz 

where the Gaussian weighted average of the frequency shift curvature is 
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and X is the half width at half maximum of the the axial density distribution 

Af(z) = c-^'/2fl" (125) 

given by X = i?V21og 2, where R = y^ksT/ muji. It is apparent from Figure 2 of [[l|] that 
we can accurately represent the gradient of the differential potential in terms of a Gaussian 
with amphtude Ag defined as 

V^z) = -AgN{z). (126) 



We then obtain Ag using (|123[),(|124|) and (|126|), to find 

erf(3i/2y2i?) 



Ag = 47ri?''ma;diff 



/^erf (3x/2i?) + {ix/R) exp (-(3x/2)2/i?2 
If we now expand (|126|) for small z we find 



47ri?^mtj^iff X 0.44. (127) 



AgN{z) ~ constant + — ( WdifiV2(27r)i/'^V(l44] ' z^, (128) 



2 



.g.,,,., , 2 

SO that the true frequency near the center of the trap is ~ 1.49i'diff- Such a correction is 
expected since the Gaussian weighting process will reduce the resulting z^diff by adding in 

negative frequencies, if the range of integration is taken beyond the turning point of the 
Gaussian iV(x). 

6.3. Results 

Figure ^ shows the result of simulating equations ( pjjp^p^l) w ith fdiff ~ 0.16 Hz, -/V(O) 



1.8 X id^^cm""^, with the relaxation time approximation ( |1Q5[ ) modelling the damping. The 
effect is very similar to that seen in the experiment In (a) and (b) the axial density profile 
of the |1) and |2) states show a clear axial segregation of the two internal states, (c) Shows the 
decay of the coherence which causes the transience of the phenomena on the timescale of the 
experiment. The initial condition is shown for comparison. Column (d) is the phase of -F(x), 
from which it is clear that the relative motion tends to smooth out the relative phase gradient 
that is caused by the differential potential Vdiff(x). 

Figure || shows the results of simulations for the experimental parameters, with the 
measured data [|l^. The timescale for decay of coherence is well described by the simple 
relaxation time approximation, with an extra factor of 2. The initiation times, relaxation times 
and amplitude of the segregation agree qualitatively with the experiment. 

Figure ^ shows the variation with Vdm, and the numerical results. We see here that 
what is described as "higher order effects" are not well modelled by the simulation for 
VdiS = 0.20 Hz. The initiation time and transience are still in good agreement but the extra 
peak in the center of iVi (x) does not emerge. 

6.4. 'Classical' Motion 

It is useful to consider the classical motion that would arise if the system could be prepared 
with identical initial state distributions as constructed by the first tt/2 pulse of the experiment. 



while also setting the coherence to zero. If coherence is neglected the initial condition (1 14) 
is merely a density perturbation from the new Boltzmann equilibrium 

n(x,p) = AAexp - cosh 



2mkBT ksT I V ksT 
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Figure 2. Time is in ms down the left column. N{0) = 1.16 X lO^^cm-^, T = 850nK, 
t'diff = 0.1. (a) and (b) are the densities of the two internal states in arbitrary units, (c) Shows 
the coherence density, with initial condition for comparison, and d) shows the relative phase 
between the two internal states = anglc(F(x)); the plot range is [0, 2tt). 



m(x, p) = exp — 



2mkBT 

The corresponding single species distributions are 



Kff(x)\ . , ^ydiff(x) 

smn 



(129) 



ii(x, p) = A/^cxp - 



P 



712 (X, P) = A/2 CXp ( - 



2mkBT 

„2 



P 



2mkBT 



Vijx) + 2miiVi(x) +m2jV2(x) 

ksT 
ksT 



(130) 

for appropriate normalizations A/i and A/2 ■ We have seen in section |6.1.1| that the interactions 
have a negligible effect on the initial condition, and similar reasoning shows the equilibrium 
segregation after the pulse is negligible for the densities used in the experiment [p. It 
follows that the classical motion towards the new equilibrium would be small, as long as 
the fluctuations caused by the perturbation are also small. 
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Figure 3. Measured N\{z), with simulations (dotted line) for Vd\s = 0.1 and peak densities 
listed at the foot of each column. Time is in ms down the left column. 

7. Spin waves and atomic motion 

In the experimental report it is asserted that the motion observed must be due to actual 
physical motion of the atoms since energy conservation would seem to prohibit spontaneous 
interconversion of the internal states [|l]|. In this situation 'motion' refers to the redistribution 
of the atomic density in the trap. The issue is whether or not it is possible to distinguish 
between two interpretations: The apparent motion is the result of either 

i) Transport of the atoms in the relative phase gradient arising from the Zeeman and mean 
field effects, or 

ii) The coherent interconversion of the internal states as a result of the interactions. 

In defence of the first picture, it has been suggested in [|l]] that the rate of segregation can be 
accounted for if the small differential potential is somehow amplified in a coherent collision 
processes which channels the thermal energy of atoms in a particular direction. In the second 
picture, the collisions between the two different internal states are interpreted as causing a 
rotation in pseudo-spin space of the spinor wave function for the two state superposition of 
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Figure 4. Measured Ni{z), with simulations (dotted line) for peak density A'^(O) = 
1.8 X 10^^ cm^'^ and u^^g listed at the foot of each column. Time is in ms down the 
left column. 



the field operators. 

It is important to note here that if the transport terms are neglected in the equations of 
motion, segregation does not occur The equations of motion for the densities are found by 
integrating the equations of motion (pj[), (p^, ( p3| ) over momentum which leads to 

7V(x) + V • (iV(x)vAr(x)) =0 (131) 

M(x) + V • (A/(x)vA/(x)) = (132) 

Fix.) + V • (F(x)vf(x)) = - iLUBi^)F{x.), (133) 

where the velocities are, for example 



/d^ppn(x,p), 
^(x) J 



(134) 



and the Ramsey frequency is ijJi?(x) = Aujc{x) + Aw^(x). 

If the transport terms are neglected the velocities vanish, and no change in the densities 
can occur — the only effect is to make a redistribution in momentum space. Inclusion of the 
transport terms then translates this redistribution in momentum space into a redistribution in 
position. In this sense, we agree with the view of that the segregation is caused by actual 
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motion of the atoms. The effect is however quite large, since it involves more forces than 
simply the classical forces induced by the gradient of the differential potential, namely the 
forces caused by the production of a spatially dependent phase of /(x, p), as seen in the 
second line of (]l5|). 

8. Comparison with other work 

While preparing this work several other papers on the JILA experiment [|l]] have appeared 
[^, ^ 1]]. The works [|], ^] use essentially the same formalism as our own, with slightly 
different degrees of approximation. In the simulations agree with experiment to the same 
quantitative degree as our own, although they obtain better results for the variation of t'dis , in 
particular we do not see higher order effects for the highest value i/jiff = 0.2 Hz as reported 
in the experiment and found by Williams et al [||]. In the simulated initiation times appear 
more rapid than seen in the experiment. Our results show better agreement for this feature, 
which may be because we have inverted the curvature averaging process used to characterize 
the differential potential. In ^ the timescale of initiation for the segregation is found 
analytically, and in [|[] it is noted that this also gives the position of the nodes of M(x). In [|^] 
a so called Landau-Lifshitz equation of motion was obtained, and solved numerically, finding 
qualitatively similar behavior, although the timescale of evolution is significantly longer than 
that observed in the experiment. 

More recent work at JILA [||] has focused on providing data for a collective mode 
analysis which has been carried out by Nikuni et al [^] using a truncated moment equation 
approach. In this experiment the differential potential is controllable in real time, so that a 
small amplitude fluctuation may be excited, and the potential gradient then set to zero for the 
remainder of the motion. The same relaxation time approximation is used to that derived here, 
and the simulations of |^] show excellent agreement with the data of Some discrepancy 
is found for the quadrupole mode in the regime where Landau damping is the dominant form 
of dissipation, and the authors note that the details of the trap are likely to be important in 
describing this process accurately. 

9. Conclusions 

We have simulated the JILA experiment using the Wigner function approach and found 
reasonable agreement with the data of We have discussed the Ramsey frequency and 
the local transition frequency for this system, and emphasized the difference of the two, 
finding that the Ramsey technique is insensitive to the decay of coherence. A relaxation time 
approximation for the experiment has been derived, including the effect of relative velocity 
which is expected to be important for higher temperature experiments. A Quantum Boltzmann 
collision term for two species with arbitrary S-wave interaction strengths was found, and 
we evaluated the effects of scattering on the coherence for non-condensed ^^Rb using a 
Boltzmann equilibrium form for the momentum distributions. It is apparent that the cubic 
terms in the collision integral depend on the phase space density and will be important for 
more degenerate regimes. 
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Appendix A. Continuum limit 

The derivation of (js^) from the terms (|2|) and (|3|) is simpHfied by the identities of the form 



ij I J, e 



ijIJ,e 



(A.l) 

so that ( |82| ) and ( ^3| ) are complex conjugates with I ^ m. The Quantum Boltzmann equation 
(Q) is then obtained by using 

7r\3 



N 



(A.2) 



so that the commutators may be evaluated and we may carry out the procedures in for 
passing to the continuum limit. This involves factorizing the operator averages and using the 
local equilibrium forms 

(Wki (x)VAnK2 (x')> = g(x - x') fhn ( , Ki ) (5ki,K2 



Ki + 5i,n (^Ki.Ks, 



(A.3) 

(A.4) 



and the continuum limit of the summation 

X (^)'AfA(K + K2-K3-K4) 



K2,K3,K4 



1 



d^K2 I d^Ks / d^Ki 



(27r)6 

X S{K + K2- K3 - K4) 



where 



1 



sin Aa 



sin Ay 



y 



sin Az 



(A.5) 



(A.6) 



and Ma(Q) is the approximate delta function arising from the factorization of operator 
averages 

Ma(Q)= J d^yj dVg(y)g(y')[g(y-y')]'c'« '''^'''' 

(This corrects an error in Eqs. (128) and (129) of [^], which led to an extra factor of (27r)^ in 
the final Uehling-Uhlenbeck collision term (132)). 

Appendix B. Integrals 

Some useful integrals for this paper are 

lQ{a,b)^ d(P d9 sine dQ Q^^^^'^'+bQcose ^g j^ 

Jo Jo Jo 
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and 



2t . 

3^— ag'+bg cos i 



where 



Defining 



Ig{a,b)= / d(t> de sine g^e"'"? (B.3) 

Jo Jo Jo 

= ^^^g(6/2V^), (B.4) 



= ^ x+ (1 + 2x2) ^^erf(x) . (35^ 



B{p) = {2TTnf [^y^^ exp (-ap^) (B.6) 
the integral 

X 5(Pi + P2 - P3 - P4)S{ei + 62 - £3 - £4) (B.7) 
is found by using the transformation Q = Pi +P2,q = Pi — P2to obtain 

/34 = ^^(^)Y d^Q J d^K J d\ J d^k cxp(-a(Q2 + q2)/2) 
X (5(q2 - k2)5(Q - K) 

/Q(a/2,0)/,(a/2,0) = W — — — . (B.8) 



'34 4 



24(27rft)3 



The integral 1234 is found in a similar fashion as 

d^Pi f d^P2 f d^P. 



1 



234 



^ J d^p^ i?(P2)i?(P3)B(P4) 



i2TTh)^ J (27rft)3 y (2t 

X <5(Pi + P2 - P3 - P4)J(ei + 62 - £3 - £4) 

V3 



47rfci3r' 



(B.9) 



Acknowledgements 

We would like to thank Rob Ballagh for his numerical expertise and Eric Cornell for informing 
us about the experiment [jlj] . Funding for this project was provided by the Marsden Fund under 
contract PVT-902. 

References 

[1] H. J. Lewandowski, D. M. Harber, D. L. Whitaker and E. A. Cornell, Phys. Rev. Lett 88, 070403 (2002). 

[2] C. W. Gardiner, A. S. Bradley, J. Phys. B 34, 4673-4687, (2001) 

[3] M. O. Oktel, L. S. Levitov, Phys. Rev. Lett. 88, 230403 (2002). 

[4] J. N. Fuchs, D. M. Gangardt, F. Laloe, Phys. Rev. Lett. 88, 230404 (2002). 

[5] J. E. WilHams, T. Nikuni, C. W. Clark, Phys. Rev. Lett. 88, 230405 (2002). 

[6] L. E. Reichl, A Modern Course in Statistical Physics (University of Texas Press, USA, 1980). 

[7] C. W. Gardiner, R Zoller, Phys. Rev. A 55, 2902 (1997). 

[8] J. M. McGuirk, H. J. Lewandowski, D. M. Harber, T. Nikuni, J. E. Williams and E. A. Cornell, Phys. Rev. Lett. 

89, 090402 (2002). 
[9] T. Nikuni, J. E. Williams and F A. Cornell. | cond-rp at/020.'i45r . 



ond-mat/Ol 



mw7 



[10] We use the raw data quoted in 
[11] D. M. Harber, H. J. Lewandowski, J. M. McGuirk and E. A. Cornell, |cond-mat/0208294 



